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Abstract — Negative imaginary (NI) systems play an important 
role in the robust control of highly resonant flexible structures. 
In this paper, a generalized NI system framework is presented. 
A new NI system definition is given, which allows for flexible 
structure systems with colocated force actuators and position 
sensors, and with free body motion. This definition extends the 
existing definitions of NI systems. Also, a necessary and sufficient 
conditions are provided for the stability of positive feedback 
control systems where the plant is NI according to the new 
definition and the controller is strictly negative imaginary (SNI). 
These general stability results captures all previous NI stability 
results which have been developed. The stability conditions in this 
paper are given purely in terms of properties of the plant and 
controller transfer function matrices, although the proofs rely 
on state space techniques. Furthermore, the stability conditions 
given are independent of the plant and controller system order. 
As an application of these results, a case study involving the 
control of a flexible robotic arm with a piezo-electric actuator 
and sensor is presented. 



I. Introduction 

Flexible structure dynamics arise in many areas such as 
flexible robot manipulators ([T], ground and aerospace vehicles 

0, atomic force microscopes (AFMs) (3], |T| and other nano- 
positioning systems |5|-[8|. In general, flexible structures can 
be modeled as infinite dimensional distributed parameter sys- 
tems f9|. However, finite dimensional approximations to these 
infinite dimensional dynamics are often used for the purpose 
of designing controllers JO), ifTO) . In designing controllers 
for these flexible systems, it is important to consider the 
effect of highly resonant modes. Such resonant modes are 
known to adversely affect the stability and performance of 
flexible structure feedback control systems [ir|-|T3l, and are 
often very sensitive to changes in environmental variables. For 
instance, a small change in the environment of the system such 
as changing temperature, can lead to significant changes in 
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the resonant frequencies of such systems. These changes in 
resonant frequencies can lead to large changes in the gain and 
phase of the system frequency response at a given frequency, 
which may lead to instability or poor performance in the 
corresponding feedback system. In addition, highly resonant 
modes lead to vibrational effects which limit the ability of 
control systems in achieving desired levels of performance in 
many applications such as precision instrumentation, optical 
systems, precision machine tools, wafer steppers, telescopes, 
and atomic force microscopes [11]. These issues arising from 
the presence of highly resonant modes in flexible structures 
motivate the need for tools to guarantee robust stability and 
performance in flexible structure control systems. 

One common solution to issues of robustness, stability, 
and performance in the control of highly resonant flexible 
structures is to use force actuators combined with colocated 
measurements of velocity, position, or acceleration llTTI - 
1131 . Here, colocation refers to the fact that the sensors 
and actuators have the same location and the same direction 
lfT4ll . Colocated control with velocity measurements, known 
as negative-velocity feedback, can be used to directly increase 
the effective damping in the system, thereby facilitating the 
design of controllers that can guarantee closed-loop stability 
in the presence of parameter variations and unmodeled plant 
dynamics lfT4ll . Similarly, a class of colocated controllers with 
position measurements, known as positive-position feedback 
controllers, where velocity sensors are replaced with position 
sensors, can also be used to increase damping in flexible 
systems as discussed in lfT2l . Also, positive-position feedback 
controllers are robust against uncertainties in resonant frequen- 
cies as well as unmodeled plant dynamics, in a similar way 
to negative-velocity feedback controllers llT2l . ifTSl . ifTSl . 

The properties of negative-velocity feedback has been stud- 
ied using passivity theory and the theory of positive real (PR) 
linear time invariant (LTI) systems; e.g., see |[T6l . ifTTI . How- 
ever, PR theory cannot be used directly when using position or 
acceleration measurements [13|. This drawback is important 
in applications to the field of nanotechnology, especially for 
nano-positioning systems, where position measurements are 
widely used; see e.g., IS, ISi-^, IM-IM- Similar issues 
also arise in application to the area of robotics where position 
measurements are also widely used. 

Lanzon and Petersen introduced a notion of negative imagi- 
nary (NI) systems in llT3l . ifTSl for the robust control of flexible 
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Fig. 1. A negative-imaginary feedback control system. If the plant transfer 
function matrix G{s) is NI and the controller transfer function matrix G{s) 
is SNI, then the positive-feedback interconnection is internally stable if and 
only if the DC gain condition, \max{G{Qi)G(Qi)) < 1, is satisfied. 



Structures with force actuators combined with position or 
acceleration sensors. In the single-input single-output (SISO) 
case, NI systems are defined by considering the properties of 
the imaginary part of the system frequency response G{juj) 
and requiring the condition j {G{juj) — G{juj)*) > for 
all UJ € (0,oo). The NI property arises in many practical 
systems. For example, such systems arise when considering 
the transfer function from a force actuator to a corresponding 
colocated position sensor (for instance, a piezoelectric sensor) 
in a Ughdy damped sti-ucture L3J, L12], US, JH], ||24J. Another 
area where the underlying system dynamics are NI, are nano- 
positioning systems; see e.g., El-llll, EH-ES. Also, the 
positive-position feedback control scheme in llT2l . 1251 . can be 
considered using the NI framework. Furthermore, other control 
methodologies in the literatures such as integral resonant 
control (IRC) [26| and resonant feedback control [53, [l28l, fit 
into the NI framework and their stability robustness properties 
can be explained by NI systems theory. 

The stability robustness of interconnected NI systems has 
been studied in lfT3l . ifTSll . In these papers, it is shown that a 
necessary and sufficient condition for the internal stability of 
a positive-feedback control system (see Fig. [U consisting of 
an NI plant with transfer function matrix G(s) and a strictly 
negative imaginary (SNI) controller with transfer function 
matrix G{s) is given by the DC gain condition 



XmaAG{0)G{0)) < 1, 



(1) 



where the notation Xmax{') denotes the maximum eigenvalue 
of a matrix with only real eigenvalues. This stability result 
has been used in a number of practical applications. For 
example in 1241 . this stability result is applied to the problem 
of decentralized control of large vehicle platoons. In [|3], 
im, the NI stability result is applied to nanopositioning in 
an atomic force microscope. A positive position feedback 
control scheme based on the NI stability result provided in 
ifTJl . ifTSl is used to design a novel compensation method for 
a coupled fuselage-rotor mode of a rotary wing unmanned 
aerial vehicle in ||29l . In ||8l, an IRC scheme based on the 
stability results provided in lfT3l . lITSI is used to design an 
active vibration control system for the mitigation of human 
induced vibrations in light-weight civil engineering structures, 
such as floors and footbridges via proof-mass actuators. An 




Fig. 2. Block diagram of a disk-drive reader head system which includes 
a flexible stmcture driven by a voice coil motor (VCM). The parameters 
R, Kt , J, and Kt, are the coil resistance, torque constant, moment of ineilia, 
and back electromotive force gain, respectively. Also, the parameters K f,C,, 
and ujn are the flexible structure gain, damping ratio, and natural frequency 
respectively. 



identification algorithm which enforces the NI constraint is 
proposed in ll30l for estimating model parameters, following 
which an Integral resonant controller is designed for damping 
vibrations in flexible structures. In addition, it is shown in ||3T| 
that the class of linear systems having NI transfer function 
matrices is closely related to the class of linear Hamiltonian 
input-output systems. Also, an extension of the NI systems 
theory to infinite-dimensional systems is presented in ll32l . 

The NI framework presented in |fT3l , IfTSl considers systems 
with poles in the open left half of the complex plane. This 
theory has been extended in |33| to include NI systems 
with poles in the closed left half of the complex plane, 
except at the origin. Also, further extensions to NI systems 
theory include the study of NI controller synthesis ll34l . ||35]| . 
connections between NI systems analysis and /i-analysis f36l, 
and conditions for robust stability analysis of mixed NI and 
bounded-real classes of uncertainties ll37l . Furthermore, the 
concept of lossless NI transfer functions is introduced in 
1 38 1, an algebraic approach to the realization of a lossless NI 
behavior is presented in ll39l . and a spectral characterization 
of NI descriptor systems is discussed in [40 1. The NI systems 
theory can be extended to nonlinear systems using the concept 
of counter-clockwise input-output dynamics as presented in 
|4r|-|43|. In |43|, a sufficient conditions under which a semi- 
linear Duhem model is counter-clockwise is given, where the 
counter-clockwise input-output system is restricted to periodic 
input signals. 

Despite generalizations of the NI systems framework pre- 
sented in ll33l . an important class of systems that cannot be 
captured by the existing NI systems framework corresponds to 
flexible systems with free body motion. These systems arise 
in areas such as rotating flexible spacecraft ll44l . rotary cranes 
[i45 l. robotics and flexible link manipulators fSSl, ll28l . P6l . 
and dual-stage hard disk drives P7l -ll50l. Flexible structures 
with free body motion lead to dynamical models including 
poles at the origin, which is not covered in earlier work on 
NI systems theory. In particular, the stability condition ([T]) is 
not well defined in the case of flexible structures with free 
body motion which results in poles at the origin, since in 
this case, the plant DC gain G(0) will be infinite. However, 
control systems involving flexible structures with free body 
motion arising in these important application areas still suffer 
from the stability and performance issues mentioned above. 
Thus we are motivated to extend the NI robust stability theory 
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developed in |fT3l . lfT31 . Il33l so that it can be applied to control 
systems involving highly resonant flexible structures with free 
body motion. 

Fig. |2] shows a block diagram of a system which includes 
a flexible structure with free body motion that arises in a 
problem of disk-drive control; see [,511. Here, a voice coil 
motor (VCM) is used to actuate the arm of the reader head. 
The free body motion of the reader head leads to a transfer 
function from the input Vi„ to the VCM output Vout which has 
poles at the origin. Furthermore, the overall system satisfies 
the NI frequency response property and includes poles at the 
origin. However, the NI stability results presented in [13], ifTSl , 
ll33l do not allow for poles at the origin and cannot be applied 
to control systems such as this disk-drive control system. 

In this paper, we present a new generalized definition of NI 
systems which allows for flexible structures with colocated 
force actuators and position sensors and with free body mo- 
tion. This definition extends the previous definitions of NI 
systems presented in ifTsl . ifTsl . ll33l . We also derive new 
generalized stability conditions for positive-feedback control 
systems involving an NI plant and an SNI controller. These 
results extend previous stability results presented in llT3l . 
lITSl . Il33l . to allow for this new definition of NI systems. 
As such, the new stability results presented in this paper 
provides definitive stability results for the positive feedback 
interconnection between an NI system and an SNI system. 

As in ifTSl . ifTSl . Il33l . the stability conditions presented 
in this paper, are given purely in terms of properties of 
the plant and controller transfer function matrices, although 
the proofs rely on state space techniques. Furthermore, the 
stability conditions given are independent of the plant and 
controller system order and can be stated without using the 
fact that the plant and the controller transfer function matrices 
are rational. 

The proofs of the NI stability results given in this paper 
and in the papers llT3l . Il33l . [[34 1 use the fact that NI systems 
can be transformed into PR systems and vice versa under 
certain technical assumptions. However, this equivalence is not 
complete as mentioned in ifTSl , lf33l . Oil- For instance, such a 
transformation applied to an SNI system always leads to a non- 
strict PR system and thus SNI systems cannot be transformed 
into strictly positive real systems. Hence, the passivity theorem 
lfT6l . ifTTl cannot be used to prove the stability of the positive 
feedback interconnection of an NI and an SNI system and 
a more sophisticated methodology is required to prove the 
stability results of this paper as well as the previous results in 
mi 1331, IMl- 

Preliminary conference versions of the stability results pre- 
sented in this paper were presented in ll52l . Il53l . However, 
in this paper, much more general versions of these stability 
results are presented in Theorem[T]and Corollaries I3J4J1 I which 
allow for the existence of free body motion in some but not 
all input-output channels. This is important since multivariable 
control systems involving flexible structures with free body 
motion usually include free body motion in some but not all 
input-output channels. This includes the case study considered 
in Section HV] which involves the control of a flexible robotic 
arm with a piezo-electric actuator and sensor along with a 



torque actuator and hub angle sensor Also, this study has not 
been considered in the previous conference versions of the 
paper. The results presented in [52] are equivalent to Corollary 
|3]and Corollary |2] in this paper, which are special cases of the 
main results reported here for the first time. 

This paper is further organized as follows: Section HI] recalls 
the existing definition for NI systems and outlines the notation 
that will be used in the rest of the paper Section Ulllintroduces 
the new generalized definition for NI systems, which allows 
for systems with free body dynamics. Also in this section, we 
present the main stability results in Theorem [T] Theorem |3] 
and Corollaries [T]|3] Section |IV] presents a case study, which 
involves a flexible robotic arm, as an application of the NI 
theory presented in this paper The paper is concluded with a 
summary and remarks on future work in Section |V] All proofs 
of the presented theorems, lemmas and corollaries are given 
in the Appendix. 

II. Preliminaries and Notation 

In this section, we recall the existing definitions of NI and 
SNI systems as given in |33l for systems with poles in the 
closed left half of the complex plane, except at the origin. 
We will use this existing definition of SNI (first introduced 
in fl5\) systems but in the next section we will present our 
new definition of generalized NI systems. We also define 
notation used to describe positive feedback interconnections 
and internal stability, which will be used to present the main 
results in this paper 

Consider the following LTI system, 

x{t) = Ax{t) + Bu{t), (2) 
y{t) ^ Cx{t) + Du{t), (3) 

where A £ R"><",S G K"^'",C G R™^", D G rx", 
and with the square transfer function matrix G(s) = C{sl — 
A)^^B + D. The transfer function matrix G(s) is said to be 
strictly proper if G(oo) = D ^ 0. We will use the notation 
' A B ] 

_____ to denote the state space realization y). 

The existing definition of NI systems states that a square 
transfer function matrix G(s) is NI if the following conditions 
are satisfied |[33l : 

1. G(s) has no pole at the origin and in Re[s] > 0. 

2. For all w > such that jco is not a pole of G(s) and 

j(GO-c.)-G(jL.)*)>0. 

3. If jujo with ojo > is a pole of G(s), it is at most a 
simple pole and the residue matrix Kq = limj^ji^g (s — 
jujo)sG{s) is positive semidefinite Hermitian. 

Definition 1: [33 1 A square transfer function matrix G(s) 
is SNI if the following conditions are satisfied: 

1) G(s) has no pole in Re[s] > 0. 

2) For all u > 0, j {G{juj) - G{juj)*) > 0. 

Now, consider a positive feedback interconnection between 
an NI system with transfer function matrix G(s) and an SNI 
system with transfer function matrix G(s) as shown in Fig. 
[T] Also, suppose that the transfer function matrix G(s) has a 
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minimal state space realization 



minimal state space realization 
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and G{s) has a 



Furthermore, it 



is assumed that the matrix / — DD is nonsingular. Then the 
closed system has a system matrix given by 



A + BD{I - DD)- 

B{I -Doy^c 



C 



BC + BD{I ~ DD)-^DC 
A + B{I - DOy^DC 



(4) 



Moreover, the positive feedback interconnection between G{s) 
and G{s) as shown in Fig. [1] and denoted [G{s),G{s)] is said 
to be internally stable if the closed-loop system matrix A in 
(|4|l is Hurwitz; e.g., see 1541 . 

III. Main results 

The main contribution of this paper is a generalization of 
the framework for NI systems presented in ll33l . Firstly, we 
introduce a new definition of NI systems that will allow for 
systems with free body dynamics. This is achieved by relaxing 
the condition requiring the absence of poles at the origin in 
the NI definition presented in |[33l . Secondly, this generalized 
definition will be used in a new set of stability conditions 
that will allow for NI systems with free body motion to be 
included into the framework of NI systems theory. Henceforth 
when a system is said to be NI, we will mean NI as defined 
below, not NI as defined in earlier papers and recalled in the 
text description above. 

Definition 2: A square transfer function matrix G{s) is NI 
if the following conditions are satisfied: 

1) G{s) has no pole in Re[s] > 0. 

2) For all w > such that juj is not a pole of G(s), 



j(G(j^)-G(j^)*) >0. 



(5) 



3) If s = juo with Wo > is a pole of G(s), then it is 
a simple pole and the residue matrix K — lim (s — 

s — >jujo 

jujo)jG{s) is Hermitian and positive semidefinite. 

4) If s = is a pole of G(s), then lim s'^G{s) = for 

S !-0 

all fc > 3 and lim s^G(s) is Hermitian and positive 

s ^0 

semidefinite. 

UnUke the NI definition presented in ll33l . Definition |2] 
allows for poles at the origin. In this case, we cannot use the 
existing stability results presented in Iil3l . ll33l . f34j . because 
the stability condition in ([T]l is not defined. The inclusion of 
poles at the origin extends the NI systems theory to include 
flexible systems with free body dynamics. In order to derive 
a new set of stability conditions that allow for systems with 
free body motion, we define the following constant matrices 
for a given to x m NI transfer function matrix G(s) : 



Go 



Gi 



Gn 



lim s^G{s), 

s — >0 



lim s 

s >Q 



lim 



G(.)- 



G2 
2 



S 

G2 
„2 



Gi 



(6) 
(7) 
(8) 



The matrices defined in (|6j-([8j carry information about 
properties of the free body motion of the system under 
consideration. These matrices will be used in stability con- 
ditions for the positive feedback interconnection of NI and 
SNI systems. Note that the DC gain condition ([T]i cannot be 
defined for an NI system with transfer function matrix G(s) 
unless G2 = Gi = 0, which reduces to the case where the 
dynamical system has no free body motion. From Condition 
4) in Definition |2] the matrix G2 is required to be Hermitian 
and positive semidefinite. Hence, it follows (e.g., see ll55l ) that 
if G2 ^ 0, it can be decomposed in the form 



G2 — J , 
where J is a full column rank matrix. 



(9) 



Now we present new conditions for the stability of a positive 
feedback control system involving an NI plant with free 
body motion. These conditions are stated using the quantities 
defined in (|6ll - ([8]). First, we define the 2m x 2m Hankel 
matrix F as 



r = 



Gi G2 
G2 



(10) 



Suppose that F ^ 0. Using the singular value decomposition 
(SVD), we can decompose the Hankel matrix F as 



F = [Hi H2] 



's 6 



























< n 


, Ui 



(11) 



j^mxn ^jjj jjjg [xiatrices U and Vi have orthogonal columns. 
Furthermore, we can decompose the fixfi matrix U1U2 using 
the SVD as 



UIU2 



u 



'Si 6 




" - T" 

Vi 







- T 
V2 _ 



(12) 



where U € M"^" and V G M"^" ai-e orthogonal matrices, 
V2 e K"^" and Si > 0. 

We now introduce some notation which will be used 
throughout the paper Given matrices X E jjmxm ^j^^ 
Y e M"^" such that det{Y'^XY) ^ 0, then the matrix valued 
function Y) is defined by 

ViX,Y) = X - XY {Y'^XYy^Y'^X. (13) 
Using this notation, we define the matrix 

Nf^V{GiO),F), (14) 
where the m x fi matrix F is given by 

F^UiV2, (15) 
and we wiU assume that det{F'^ G{0)F) ^ 0. 

We will use the following condition in the theorem which 

F^G(0)F < 0. (16) 



follows; 
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Also, for the case in which Nf is positive semidefinite, we 
will use the condition 

/ - NjGoNj - NjGiJiJ^jy^J^G'^Nj > 0. (17) 

Moreover, for the case in which Nf is negative semidefinite, 
we will use the condition 

det(/ + NfGoNf + NfGiJ{J^ jy^J^G^ Nf) ^ 0. (18) 

Here, Nf — {—Nf)^ and matrices Gi,Ga, J, Nf and F are 
defined in (|7]l, Q, ( fT4l i. and (fTsT i respectively. Also, 
denotes the square root of a positive semidefinite matrix. 

The following theorem is our main stability result for the 
case in which G2 7^ 0. This corresponds to the case of plants 
which have free body motion without friction being present. 

Theorem 1: Suppose that the square transfer function ma- 
trix G{s) is strictly proper and NI with G2 ^ 0, and the 
transfer function matrix G(s) is SNI. Also, suppose that the 
matrix F^G{0)F is non-singular. If Nf is positive semidef- 
inite, then the closed-loop positive-feedback interconnection 
between G{s) and G{s) as shown in Fig.[T]is internally stable 
if and only if conditions ( fT6] l and ( fTTj i are satisfied. Further- 
more, if Nf is negative semidefinite, then the closed-loop 
positive-feedback interconnection between G(s) and G{s) is 
internally stable if and only if conditions (fTST i and (fTSl l are 
satisfied. 

The proof of this theorem is presented in the appendix. 

Theorem [T] gives a stability result for the positive-feedback 
interconnection of an NI system including free body motion 
and an SNI system for the case in which G2 ^ 0. We now 
present a corollary to this theorem which considers the special 
case in which none of the free body modes of the plant have 
frictional force present; i.e., Gi — 0. In order to present this 
corollary, we define the matrix N2 as follows: 



N2=r{GiO),J), 



(19) 



where we assume that the matrix J^G(O) J is non-singular 

We will use the following condition in the next corollary, 
which corresponds to condition (fT&t in Theorem [T] 



J^G(O) J < 0. 



(20) 



Also, for the case in which N2 is positive semidefinite, we will 
use the following condition which corresponds to condition 
(fTTI l in Theorem [T] 



I-NiGoNi > 0. 



(21) 



Moreover, for the case in which A^2 is negative semidefinite, 
we will use the following condition which corresponds to 
condition (fTSl l in Theorem [T] 



dct(/ + N2G0N2) ^ 0, 



(22) 



where N2 = (— -/V2)^- 

Corollary 1: Suppose that the transfer function matrix G(s) 
is SNI and the strictly proper transfer function matrix G(s) 
is NI with Gi — {) and G2 ^ 0. Also, suppose that the 
matrix j'^G{0)J is non-singular If N2 is positive semidef- 
inite, then the closed-loop positive-feedback interconnection 
between G(s) and G(s) is internally stable if and only if 



conditions (|20] | and (I2TI 1 are satisfied. Furthermore, if N2 is 
negative semidefinite, then the closed-loop positive-feedback 
interconnection between G(s) and G(s) is internally stable if 
and only if conditions (|20] | and (l22T i are satisfied. 

The following theorem imposes some extra conditions on 
the matrix G2 which enables us to relax the sign definiteness 
condition on the matrix A^2- This then leads to a simplified 
stability condition. 

Theorem 2: Suppose that the transfer function matrix G(s) 
is SNI and the strictly proper transfer function matrix G(s) 
is NI with Gi — and G2 7^ 0. Also, suppose that 
7V(G2) C A/'(Go ), where 7V(-) denotes the null space of a 
matrix. Then the closed-loop positive-feedback interconnec- 
tion between G(s) and G(s) is internally stable if and only if 
condition (|20] | is satisfied. 

An important special case of the above theorem is the 
following corollary. 

Corollary 2: Suppose that the transfer function matrix G(s) 
is SNI and the strictly proper transfer function matrix G(s) is 
NI with Gi = and G2 > 0. Then, the closed-loop positive- 
feedback interconnection between G(s) and G(s) is internally 
stable if and only if G(0) < 0. 

In the SISO case, as in the example shown in Fig. |2l the 
effect of free body motion may lead to a double integrator in 
the plant model. However, if frictional force is present, then the 
effect of the free body motion leads to only a single integrator. 
In the general MIMO case, this corresponds to the situation 

in which G2 = lim s^G(s) = 0. In Theorem |3] Theorem |4] 

s — >-o 

and Corollary |3] we consider cases which corresponds to free 
body motion with frictional force present. As in Theorem [T] 
this case allows for fact that the free body motion may not be 
present in all input-output channels. 

In order to present Theorem, |3] Theorem S] and Corollary 
suppose that Gi 7^ 0. Then we consider the following SVD 
decomposition of the matrix Gi defined in (|7]i: 



Gi — [Fi F2] 



'S2 


0' 

















(23) 



where ^2 > 0, and the matrices Fi — F1S2 and Vi have 
orthogonal columns. Also, we define the matrix Ni as follows: 



Ni^V{G{0),F, 



(24) 



where the matrix F^G{0)Fi is assumed to be non-singular. 

We will use the following condition in Theorem |3] and 
Corollary |3] which corresponds to condition (IT6^ in Theorem 

m 



FlG{0)Fi < 0. 



(25) 



For the case in which A^i is positive semidefinite, we also will 
use the following condition which corresponds to condition 
([TtI i in Theorem [T] 



I-N^GqN^ > 0. 



(26) 



Moreover, for the case in which Ni is negative semidefinite, 
we will use the following condition which corresponds to 
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condition ( fTST i in Theorem [T] 



dct(/ + iViGoA^i) ^0, 



(27) 



where iVi = (-iVi)^. 

Theorem 3: Suppose that the transfer function matrix G{s) 
is SNI and the strictly proper transfer function matrix G{s) 
is NI, with G2 — and Gi 7^ 0. Also, suppose that the 
matrix F^G{0)Fi non-singular If A^i is positive semidef- 
inite, then the closed-loop positive-feedback interconnection 
between G{s) and G{s) is internally stable if and only if 
conditions (|25] l and ( |26] l are satisfied. Furthermore, if A^i is 
negative semidefinite, then the closed-loop positive-feedback 
interconnection between G{s) and G{s) is internally stable if 
and only if conditions (IZST i and dZTl i are satisfied. 

The following theorem imposes some extra conditions on 
the matrix Gi which enables us to relax the sign definiteness 
condition on the matrix A^i. This then leads to a simplified 
stability condition. 

Theorem 4: Suppose that the transfer function matrix G{s) 
is SNI and the strictly proper transfer function matrix G{s) 
is NI with G2 = and Gi 7^ 0. Also, suppose that 
■^{GJ) ^ Af{GQ). Then the closed-loop positive-feedback 
interconnection between G(s) and G(s) is internally stable if 
and only if condition (|25] l is satisfied. 

The following corollary presents an important special case 
of the above theorem. 

Corollary 3: Suppose that the transfer function matrix G(s) 
is SNI and the strictly proper transfer function matrix G(s) is 
NI with G2 = and Gi invertible. Then, the closed-loop 
positive-feedback interconnection between G(s) and G(s) is 
internally stable if and only if G(0) < 0. 

Remark 1: In the case where G2 = and Gi = 0. This 
case is equivalent to the existing stability results presented in 
lfT3l . II33I . Il34l . In this case, the stability condition reduces to 
Amax(G(0)G(0)) < 1. This condition can be obtained from 
(l26T l using the fact iVi — G(0) in this case. Also, we require 
the assumption G(0) > 0. Hence, 

I-N^GoN^ > 0, 
«=>7Vf ^ - Go > 0, 
max (G(O)Go) < 1. 

Note that using a similar argument to the proof of Theorem 
[3] we can obtain a similar result under the assumption that 
G(0) < 0. 

IV. Case Study: Control of Flexible robotic arm 

In this section, we present an application of the stability 
results presented in this paper to the control of a flexible 
robotic arm system. The robotic arm is pinned to a motor 
at one end with a piezoelectric laminate film attached to the 
arm. For the purposes of modeling the flexible robotic arm, 
we use an equivalent slewing beam model as depicted in Fig. 
m see ||56l . Here, L represents the length of the robotic arm, 
xi is the distance from the center of the motor to the start 
of the piezoelectric film, and X2 is the distance from the 
center of the motor to the far end of the piezoelectric film. 
The motor allows the robotic arm to traverse in the vertical 



plane. Two piezoelectric patches are attached to the arm on 
either side. Here, one piezoelectric patch acts as an actuator 
while the other is a sensor The robotic arm system has of two 
inputs and two outputs: the inputs are the voltage Va applied 
to the piezoelectric actuator and the torque r applied by the 
motor, whereas the outputs are the voltage Vs produced by the 
piezoelectric sensor and the motor hub angle 9. The fact that 
this system involves colocated "force" actuators and "position" 
sensors indicates that the system will be NI; e.g., see l,13J . 



A. Mathematical model for the robotic arm 

The beam in Fig. [5] is modeled using the Bernoulli-Euler 
equations of motion for a beam with actuating and sensing 
piezoelectric elements as in |56|: 



EI^y(x,t) 



CaVa{x,t) 



(28) 



Here, E is Young's modulus and / is the second moment 
of inertia of the beam, p is the density of the beam, and 
A is the area of the composite beam. If the thickness of 
the piezoelectric films are comparable to the thickness of the 
beam, then the products EI and pA would be different in 
the laminated and non-laminated areas of the beam. However, 
since piezoelectric films used in practical applications are often 
thin compared to the thickness of the beam, these differences 
will be neglected. Assuming that the products EI and pA are 
uniform over the length of the beam simplifies the modeling 
procedure. 

Now we consider various boundary conditions in modeling 
the beam. These are given as 



2/(0, t) = 0, 
52 



93 



(29) 

y(0,t) + r(t) =0, (30) 
(31) 
(32) 



Here, Mt and It are the mass and inertia of the tip, which will 
be neglected in this paperAlso, ( |29] | represents the inability 
of the motor joint to undergo transverse motion. As in ll56l . 
the time domain beam equation ( l28l ) with boundary condi- 
tions (|29]l-(l32]i can be transformed into an equivalent Laplace 
domain representation as 

Y"'\x, s) - /3^Y{x, s) = 9^^Y^M. (33) 




Fig. 3. Schematic of a slewing beam equivalent of a robotic arm. 
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with boundary conditions 

Y{0,s)=0, 

EIY"{0, s) - IhS^Y'iO, s) + t{s) 
EIY"{L, s) + Its^Y'{L, s) = 0, 
EIY"'{L, s) - Mts^Y{L, s) = 0, 

where the primes indicate spatial derivatives and 

pAs^ 



computed as 



0, 



(34) 
(35) 
(36) 
(37) 



EI 



Note that ( l33T l is the Laplace domain equivalent of 
theBemoulli-Euler beam equation with V^{-) as a forcing 
input. Together, (l33Tl-(l37]i represent a set of linear ordinary 
differential equations with mixed boundary conditions: two at 
X = Q and two ?A x — L. A state space representation for the 
system can be formed from equations (I33t-(l37t as in Ii56l : 



-Y'{x,s)- 




-0 


1 





0- 




- Y{x,s) - 


Y"{x,s) 










1 







Y'{x,s) 


Y"'{x,s) 













1 




Y"{x,s) 


Y""{x,s)_ 

















Y"'{x,s) 



+ 



CaVajs) 

EI 



^5'(a:-x0(-l)' 



(39) 



where S{-) represents the Dirac delta function. The equation 
( [39] l can be written in the general form 

Z'{x,s) = AZ{x,s) + BU{x,s), (40) 
the solution to which is given by 

Z(x,s) 

= e^'^ZiO, s) + [^e-*(^~^i'B - Ac-^'^-^^'bI ^lii^iif) . (41) 
L J EI 

Once the boundary conditions Z{Q, s) and Z{L, s) are known, 
dJU will depend upon three conditions for x, namely, < a; < 
xi, xi < X < X2, and X2 < x < L. For further details see 



B. Infinite Dimensional Transfer function Model 



Here, we present the input-output relationship between the 
two inputs Va and t, and the corresponding collocated outputs 
Vs and 6 in the form of the transfer function matrix. 



0{s) 
Vs{s) 



Gis) 



r{s) 



(42) 



where G{s) 



and each of the ele- 



Gr,vAs) Gv^yAs) 
ments of this transfer function matrix is an infinite dimensional 
transfer function defined in terms of transcendental functions 
of /3. Indeed, each of the four transfer functions in (l42t can 
be written as a ratio of numerator and denominator functions 



Gr,e(s) = 

Gr.vAs) = 



N^As) ^ Y'{G,s) \ 
D(s) T A(s)=o 

D{s) Va It(s) = o' 

N^y^ {s) _ CAY'{x2,s)-Y'{x^,s))\ 



(38) Gv,,vM 



D{s) t{s) IVa(3)=0 

Nv^,Vs{s) _ Cs{Y'(x2,s)-Y'{x^,s)) 



(43) 
(44) 

^ 1 

(45) 



D[s) 



Va{s) 



lT(s)=0 

(46) 



Here, 



D{s) ^ApEI{pA{cos{l3l)smh{l3l) - cosh(/3/) sm(/?0)) 

- 4:AEIIh(l + cos{A)) cosh(/3Z)), (47) 

where //; is the hub inertia. Also, the functions 
Nrfi{s),Nv^,e, Nr,v^ (s), Ny^y^ (s) are given by very compli- 
cated expressions which can be found in equations (26)-(28) 
in |56|. 

We now compute the transfer functions in (l43ll-(l46ll for the 
case where the piezoelectric actuators and sensors span the en- 
tire length of the beam. This corresponds to the substitutions: 
xi = and X2 = L. The resulting transfer functions have 
been verified in f57l for an experimented robotic arm system. 

Despite the fact that we have not defined the NI property 
for infinite dimensional transfer functions, we will provide 
some calculations which indicate that the infinite dimensional 
transfer function matrix G{s) defined in (|42] |-(|46]| satisfies 
the NI conditions given in Definition ID Since the infinite 
dimensional transfer function matrix G{s) is actually a tran- 
scendental function of /3(s), G{s) — G(/3(s)), then Condition 

2) in Definition [2] is equivalent to the condition 

j-(G(/3(jL.))-G(/30w))*) >0 (48) 

for all w > where /3(s) is given by (|38] |. Indeed, it is 
straightforward to verify from the formulas for the transfer 
function matrix that j{G(AA)) - G{AA))*) 

for all w > 0. Also, the function D{s) given in (|47] | has an 
infinite numbers of roots. However, we can check Condition 

3) in Definition 2 for a finite number of these roots on 
the imaginary axis. To do so, we have calculated the first 
eleven joj-axis roots of D{s) numerically. At each of these 
roots s = jwo, the corresponding residue matrix K{j uiq) — 

lim (s — jujo)jG{s) is calculated using L'Hopital's rule as 

s — >jua 

follows: 



lim (s - juJo)jG{s) 



= lim (s - jujo)j 

S >JUJO 



D(s) 
D{s) 



Ny^.ejs) 
Dis) 

D{s) 



■ J 



D'Uujo) D'(juo) 



(49) 



where D'{juj) denotes the first derivative of D{jA with 
respect to uj. 

In this case study, the parameter values for the robotic arm 
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are taken to be the same as used in ||57l . These parameter 
values are shown in the Table |T] 

TABLE I 

Robotic arm parameter values. 



Parameter 


Value 


Unit 


Hub inertia, Ih 


0.0348 


N - s '-' 


Beam length, / 


2 


m 


Volumetric mass density, p 


2712.6 


Kg/m 


Cross sectional area, A 


483.87x10"" 




Young's Modulus, E 


69.0x10" 


N/m-' 


Area moment of inertia, I 


1.63 X 10"'' 




Coupling Coefficient fcai 


-0.340 




Capacitance, C 


68.35 


p.F/m-' 


Thickness ts 


3.05x10"* 


m 



180p 
160- 
= 140- 




100 200 300 400 500 

^ (rad/scc) 

Fig. 4. Plot of the log of the minimum eigenvalue of the matrix 

D' {ju))'^ K{jui) + 7D(jrcj)^ versus frequency u). 



Table HI] shows the calculated roots of D{s) and the min- 
imum eigenvalue of the corresponding residue matrix given 
in ( |49] l. Also, the matrix G2 



0.14 




lim s G(s) is found to be 

: VO 



which is positive semidefinite. These results show 



TABLE II 

The minimum eigenvalues of the residue matrix corresponding 
to the first ten resonant modes of the infinite dimensional 
plant transfer function matrix. 



n 


Root s = ju) (rad/s) 


Minimum eigenvalue 











1 


3.395326441 


0.1434 


2 


9.501801884 


0.2553 


3 


17.08210071 


0.1320 


4 


29.32863976 


0.0361 


5 


47.01240951 


0.0142 


6 


96.84550724 


0.0049 


7 


128.7332003 


0.0034 


8 


165.2195349 


0.0025 


9 


206.2898971 


0.0019 


10 


251.9420283 


0.0015 



that the infinite dimensional transfer function matrix G{s) 
satisfies the conditions of Definition |2] at least for the first 
ten resonant modes. 

Also in the Fig. |4] we plot the log of the minimum 
eigenvalue of the matrix D'{ju!)^K{juj) + 'yD{juj)'^ as a 
function of frequency ut, where 7 is a positive constant. This 
plot also indicates that the residue matrix defined in ( |49t 
will be positive semidefinite at the system poles within the 
frequency range of interest. 



C. Approximate Finite-dimensional Transfer Function Matrix 

The transfer functions in (|43]i-(|46]l are irrational functions 
of s. We now approximate these transfer functions by rational 
functions in s in order to design a suitable controller for the 
robotic arm system and to simulate its performance. Various 
methods such as the Maclaurin series expansion presented 
in |58|, the Rayleigh-Ritz method |59|, and the assumed 
modes method |59| are available in literature for the finite 
dimensional approximation of infinite dimensional systems. 
Here, we adopt a partial fraction approach to obtain a finite 
dimensional approximation of G{s). This method is similar 



to the assumed modes technique described in |,59J . The finite 
dimensional model can be written as 



Gf{s) = 



GfT,e{s) GfVa,eis) 

Gfr.Vsis) GfVa,vAs) 

ra ^ r b; 

^ k 



+Pi 



where D{s) in (|43]|-(|46]| is approximated by 

n 



(50) 



(51) 



i=0 



Here, jpo...jpn are the first n ju-axis roots of D{s) and the 
coefficient matrices C,; = , 
partial fraction expansion method. 



The constant k is chosen so that 



are computed using the 
That is, 

Nr,e{jPi) Nv^,o{jPi)' 

NT,V,{jPr) Nva,vAjP^) 



(52) 



^(j^o)=fcn(- 



(53) 



where ujq is such that Juq is not a root of D{s). We consider 
the first four resonant modes; i.e., n — 4. The corresponding 
coefficient matrices were computed using (|52] | and were found 
to be 



Go 



G2 



C3 



C4 



;Ci 



0.14 0" 
0_ 

12.1185 
-0.2165 X lO"'' 

11.0959 
11.724 X 10"^ 

3.0791 
-1.4638 X 10" 



3.0907 3.5573 x lO""* 
3.5573 X lO"'' 2.3500 

-0.2165 X 10"" 
44.1415 
11.724 X 10" 
911.9085 

-1.4638 X 10"* 
17438.04 



(54) 



and the constant k was calculated using (|53] l to be fc = 
6.6667 X 10"*. Also, the poles were computed to be pq = 
0,pi = 3.4; p2 = 9.5; P3 = 17.08; p4 = 29.33, which are the 
first four resonant modes given in Table |II] 

It is clear that the finite dimensional model Gf{s) in 



dlOll is NI, since j{G{juj) - G{juj)*) = for all a; > 
and the coefficient matrices Co, Ci, C2, C3, C4 are positive 
semidefinite which implies that Condition 3) in Definition |2] is 
satisfied. Also, G2 = lim s^G'(s) = Co > 0, which implies 

s >0 

that the Condition 4) in Definition |2l is satisfied. 



D. Controller design 

According to Theorem [1] if a plant is NI, any SNI controller 
which satisfies the conditions of Theorem [T] will stabilize the 
system. The fact that the robotic arm plant involves colocated 
"force" actuators and "position" sensors indicates that this 
plant should be NI. In particular, the finite dimensional approx- 
imation to the robotic arm model derived in Subsection llV-CI 
was shown to be NI. We will now use this finite dimensional 
model to design a controller for the system. First, we compute 
the matrices G2, Gi, and Go in ©-([Sl, for the finite dimension 
approximate system (ISOl l to obtain 



G2 = 



Go 



>0; 



Gi = 



(55) 



'0.14 0' 



"0.41253083 0.0000319 ' 
0.0000319 0.15672805 

This implies that we can use Corollary [T] to guarantee the 
stability of the positive feedback interconnection between the 
plant and an SNI controller 

In this case study, an integral resonant controller (IRC) is 
chosen to stabilize the system; e.g., see [13|. An IRC is a first 
order controller which takes the form 



G(s) = (s/-r$)-ir- A. 



(56) 



This controller is SNI if F > 0, $ > and A is a symmetric 
matrix 1 13|. Now, we chose the controller matrices F > 0, $ > 
and A such that the conditions of Corollary [T] are satisfied. 
We choose the controller matrices as follows: 





"35 15" 




r = 






$ = 




15 20 






"4.2900 


■ 




A = 













2.22 





0.745 
0.521 



0.521 
1.021 



(57) 



This leads to a controller DC gain matrix of G(0) = 
-2.2029 -1.0650l 

. To check the stability conditions in 
-1.0650 -0.697lJ ^ 

Corollary [T] we first compute the matrix J in ^ using 

^ '0.3751" 
G2 in (|55ll. This yields J 







N2 in (fT9] l will be given as 



Also, the matrix 





-0.182252 



which is 



negative semidefinite. Then we conclude det(J + A'^2GoiV2) 
ri.000000025 0.000000003^ 



det 



/ 0, where iV2 = {-N2 



1/2 

5.390603 _ ^ 
0.309908135 < 0. Thus, the conditions of 



0.000000003 
Also, J'^G{0)J = 
Corollary [T] are satisfied. 

To verify the performance of the closed loop system, we 
plot the step response of this system corresponding to a step 
change in the reference position of the robotic arm; see Fig. |5] 
This step response is shown in Fig. |6] Also, the corresponding 
response of the piezo sensor output Vs is shown in Fig.|7] Here, 
the step responses were calculated using finite dimensional 




Fig. 5. Block diagram corresponding to a step change in the robotic arm 
reference position. 




— Four Modes 

- - Three Modes 

- -Two Modes 

One Mode 



100 120 

Time (sec) 



Fig. 6. Position response of the robotic arm system corresponding to a unit 
step change in reference position. Here, different numbers of modes are used 
in the plant model. 



plant models defined in (|50]|-(|54]| for different numbers of 
modes. 

Note that the controller parameters in (l57i were chosen 
by process of trial and error to obtain good closed loop 
performance. An alternative approach, which would be useful 
in the case of a more complicated SNI controller structure, 
would be to use an optimization procedure to obtain the 
controller parameters; e.g., see 11601 . 



0.4 
0.3 

0) 0.2 

T3 

E 

< 
-0.1 
-0.2 

-0.3L- 



— Four Modes 

- - Three Modes 

- - -Two Modes 

one Mode 



80 100 120 

Time (sec) 



Fig. 7. Piezoelectric response sensor output Vs response con'esponding to 
a unit step cliange in reference position. Here, different numbers of modes 
are used in the plant model. 
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V. Conclusion 

In this paper, new stability results for the positive-feedback 
interconnection of negative imaginary systems have been de- 
rived. A new NI definition is presented, and a generalized NI 
lemma is derived. The new NI definition allows for systems 
having free body dynamics to be considered as NI systems. 
This work can be used in controller design to allow for a 
broader class of NI systems than considered in previous work. 
The application of the main results in this paper has been 
illustrated via a case study involving the control of a flexible 
robotic arm. 
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VI. Appendix 

In this appendix, we prove our main results using a partic- 
ular state space representation of the plant transfer function 
matrix G{s). The first stability result will be presented in 
Theorem |5] We will then use Theorem |5] to prove the result 
presented in Theorem [T] We also present a number of corol- 
laries which will be used to prove remaining results of the 
paper. 

We first consider an NI square transfer function matrix G(s) 
with a minimal state space realization of the form 



x{t) = Ax{t)+Bu{t), 
V{t) = Cx{t), 

where, A e M"^", B e R"^™, C G ]R™^'\ 



and 





'A, 





" 




'Bi' 


A = 





A2 





; B = 


B2 










As_ 







C^[Ci 



(58) 
(59) 



C2 C3] , 
(60) 



'-,A3 e M2fex2fe^^^ g 

A\ is nonsingular, 
(61) 



[C'3a ^3;,] , C3a G W^^^ ^C-^h £ 

A2 0, and 



A. 



/fcxfc 





Remark 2: We allow any of the matrices in the model 
to have zero dimensions. In sequel, any matrix a with zero 
dimension is regarded as being of full rank. 

The corresponding transfer function matrix for the state 
space realization (I58]l-(l60ll is given as follows: 

G{s) =Ci(s/ - Ai)-^B^ + C2{sl - A2)-^B2 
+ C3isI -Asy'Bs 

n I T A \-^U < C2B2+C3B3 CsaBsb 

=6i(s/-Ai) Bi H \ 5 — . 

s 

(62) 

The following lemma gives expressions for the quantities 
Go , Gi and G2 in ©-(HJ in terms of the state space realization 

dSSi-dSoli. 

Lemma 1: Suppose that G(s) has a minimal state space 
realization (|58]|-(|60]|. Then the quantities Gq,Gi and G2 
defined in ©-(lUl are given as follows: 



G2 = lim s^G{s) = C3aB3b, 

S ^0 



Gi = lim s(G(s) 

s — >0 

Go = lim (G(s) - 

s — >0 



G2 

G2 



G2B2 + G3B3, 



Gi 



) = -C,A-'B,. 



(63) 
(64) 
(65) 



Proof: This lemma follows immediately from 
The following lemmas give some useful properties of the 



minimal state space realization (|58]|-(|60]|. 

Lemma 2: Suppose that the transfer function matrix G(s) 
has a minimal state space realization (l58]l-(l60ll. Then, the 
matrix [G2 G3a] is of full column rank, and the matrix 

is of full row rank. Also, the subsystem with realization 



B: 



36 



Ai 


Bi 


. ^1 






is minimal. Furthermore, m > k and to > 712. 



Proof: Since the state space realization 



A 



C 



B 







IS 



minimal, the pair {A, G) is observable and the pair {A, B) is 
controllable. Also, the corresponding observability matrix is 
given by 





C 




Gi 


G2 


G3a 


G36 




GA 




GiAi 








G3a 


0{A,G) = 


CA^ 




GiAl 

















CiA'l-^ 












(66) 

Since the pair (A, G) is observable, it follows that the observ- 
ability matrix 0{A, G) in (|66] | is of full rank. This implies that 
the pair [Ai^Ci) is observable. Also, since the observability 
matrix 0{A, G) is of full rank, it follows that G2, G3Q and 
[G2 G3q] are of full rank. Furthermore, we must have m> k 
and m > n2- Similarly, since the pair {A,B) is controllable, 
it follows that the corresponding controllability matrix 

C{A,B)^[B AB A^B ■■■ A'^-^B] 



"Bi 


AiBi 


AjB, • 


• A'^^^Bi- 


B2 











Bsa 


Bsb 








.Bzb 












is of full rank. Hence, the pair {Ai,Bi) is controllable and 

'B2' 



the matrices B2, B31, and 



Bsb 



are of full rank. Also, 



since the pair (Ai.Gi) is observable and the pair {Ai,Bi) 



is controllable, it follows that 



Ai 


Bi ' 


. ^1 






is a minimal 



realization. 

Lemma 3: Suppose that the transfer function matrix G(s), 
with the minimal state space realization (l58]l-(l60ll. is NI. Then, 
there exist symmetric matrices Pi > 0, P2 > 0, and matrices 
Li, W such that 

PiAi + AlPi = -LlLi 
PiBi - AlG^ = 



P2S3& = G 



3a I 



(67) 
(68) 
(69) 



GiBi + bJ^GJ + G2B2 



G3i33 + bTgT 



= w'^w. 



Furthermore, 



(70) 
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Pi 



and 



— C2B2 + -BJcJ + C^B^ 
= [W'^ + CiP{^Ll) {W- 



(71) 



LiPi-'C^) > 0. 

(72) 



Proof: Consider the transfer function matrix G{s) with 
the minimal state space realization (|58]l-(l60ll. Also, define the 
transfer function matrix R{s) — sG{s). Using ( l62T i. it follows 
that 



R{s) ^sCi{sI ^ Ai)-^Bi 



3b 



= CiAi{sI 
+ C2B2 - 

This implies that 









Dr 



CxAi 



~A,)-'B^ 
- C3B3. 
R{s) has 



CiB, 



a state 



where Ar = 

Csa"] and Dr = C1B1+C2B2+C3B3. Using Lemma 



'Ai 6 




' Bi 




, Br = 







Bsb 



- C2B2 + C3B3 

S 

(73) 

space realization 
' 1, 1 



[2] it follows that the state space realization 









Dr 



IS 



minimal. 

We now show that R{s) is positive real; e.g., see page 47 in 
ifTTl for a definition of positive real transfer function matrices. 
Since G{s) is NI, it follows that j {G{juj) - G{juj)*) > 0, for 
all > such that ju is not a pole of G{s). Then given any 
such uj > 0, Rjjuj) + R jjuj)* ^ juj {G{juj) - G{juj)*) > 
0, and {R{jio) + R{juj)*) > 0. This implies that R{-juj) + 
R{-juj)* > for all a; > 0, which implies that R{juj) + 
> for all cl> < such that juj is not a pole of G{s). 
Hence, R{juj) + -R(jw)* > for all uj e (-00, 00) such that 
ju is not a pole of G(s). 

Now, consider the case where jujQ is a pole of G(s) and 
Wo = 0. In the case where GsaSsb = 0, the transfer function 
matrix R{s) = Gi Ai(s/ - + GiBi + Ga^a + G3-B3 

will have no pole at the origin. This implies that i?(0) is finite. 
Since R{juj) + R{juj)* > for all w > such that juj is not 
a pole of G(s) and R{jio) is continuous at w = 0, this implies 
that i?(0) + i?(0)* > 0. In the case where CsaBab 7^ 0, the 
transfer function matrix R{s) is as given in ( f73T l. Since G(s) is 
NI, then lim s^G(s) > which impHes that lim sR{s) > 0. 

s ^0 s >0 

If jujQ is a pole of G(s) and luq > 0, then G(s) can be 
factored as 



■F{s), which according to the definition for 
NI systems implies that the residue matrix Kq — a^^O'^o) is 
positive semidefinite Hermitian. Hence, F(jwo) = F^jujo)* > 
0. Now, the residue matrix of R{s) at j'wq with wq > is 
given by. 



lim (s - juja)R{s) = 



lim (s - jwo)sG(s), 

S >jhJo 

lim {s ~ juja)s^-]- — 7 



Fis), 



which is positive semidefinite Hermitian. Hence, we can 
conclude that R{s) is positive real; see page 47 in |17|. 

Using the KYP lemma (e.g., see Lemma 3.1 in 1 17]), it now 
follows that there exist matrices Pr > 0,L and W such that 



PrAr ~\~ A^ Pr — L 

PrBr - GJ = - L^W, 



Dr + DJ =W^W. 



(74) 
(75) 
(76) 



If we write Pr = 
from (l75]l that 

>i 



Pi 
^12 



P12 
P2 



and L = [Li L2] , it follows 



^12 



P12 Ai 
P2 

LI 

PiA, + AjPi AlP^2 
PI2M 



Al 




Pi 

P^ 
^12 



P12 
P2 



[Li L2] 



L2 Li 



LJL2 



L2 L2 



(77) 



>1 


0' 




' Bi 




ClAt 









P2 




B^b 




^3a 








Hence L2 — and since Ai is a nonsingular matrix, it also 
follows that P12 = 0. Also, dTTl ) impUes that Wl\ is satisfied. 
From (|75] l, it follows that 



which implies ( |68] | and ( |69] l. Lemma |2] implies that Baf, is of 
full rank and hence, ( |69] l implies that (fTTT i is also satisfied. 

From ( |76] l, it follows that (fTOl i holds. Also, using (|68] l, we 
can write _Bi as, 

Bi = Pf ^(Af G?" - LlW). (78) 
Substituting (|67]l and (|78ll into dTOll, it follows that 

C2B2 + -bJgJ + G3B3 + bJc^ 

= (VK^ + GiPi^^Lf ) (VF + LiP^^Cl) > 0. 

Using ( |64] | in Lemma [T] this implies ( |72] |. This completes the 
proof. ■ 
We now consider an SNI transfer function matrix G(s) with 
a minimal state space realization 

±{t) = Ax{t) + Bu{t), (79) 
y{t) ^ Cx(t) + Du{t), (80) 

where A e R"^«, B G R"^", G G R"^", and D e M™^". 

The following theorem provides a necessary and sufficient 
condition for the stability of the positive-feedback intercon- 
nection between the NI transfer function matrix G(s), with 
state space realization (|58]|-(|60]|. and the SNI transfer function 
matrix G(s), with state space realization (l79]l- (l80l l. In order 
to present this theorem, we define the following matrix 



7V-P(G(0),[G2 Gsa]) 



where 



CI 

^3a 



G(0)[G2 G3a] 



(81) 



(82) 
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is assumed to be non-singular. 

We will use the following condition in the theorem which 
follows: 



G(0) [C2 C3a] < 0. 



(83) 



Also, for the case in which N is positive semidefinite, we will 
use the condition 

/ + N^CiA^^BiN^ - N^CsbPi^Cj^N^ > 0, (84) 

where P2 — C^aP3bi^3bB^b)^^ ' which is symmetric and 
positive definite according to Lemma[3] Moreover, for the case 
in which N is negative semidefinite, we will use the condition 

det(/ - NCiA^^BiN + NCabP^^^C^^N) ^ 0, (85) 

where iV = . 

Theorem 5: Suppose that fc 7^ and the matrix S in 
(l82T l is non-singular Also, suppose that the transfer function 
matrix G{s)^ with the minimal state space realization (l58l l- 
(|59] |, is NI and the transfer function matrix G{s)^ with the 
minimal state space realization (l79]l-(l80ll. is SNI. If iV is 
positive semidefinite, then the closed-loop positive-feedback 
interconnection between G(s) and G(s) is internally stable if 
and only if conditions ( l83T l and ( |84] | are satisfied. Also, if N is 
negative semidefinite, then the closed-loop positive-feedback 
interconnection between G(s) and G(s) is internally stable if 
and only if conditions (|83] | and dSSl l are satisfied. 

We now present some lemmas which will be used in the 
proof of Theorem |5] 

The following lemma follows directly from the proof of a 
result presented in ll33l . 

Lemma 4: (See the proof of Lemma 6 in fSS^) Suppose that 
the transfer function matrix (^(s), with minimal state space 
realization (|29ll-(|80ll, is SNL Then, D = 5^, det(A) ^ and 
there exists a matrix P — P^ > such that 



AP- 



P-^A^ < and B 



AP-'C 



(86) 



Lemma 5: Suppose that the transfer function matrix G(s) 
with the minimal state space realization (l58]l-(l60ll is NI. Then, 

there exists an invertible matrix Rd such that 



Pd 



Also, if X e A/" 



Bsb 

|, then X e M{B3a 



P2 G|J,), where the matrix P2 is defined as in Lemma |3] 



. It follows that 



G,^ 



Proof: Suppose that a; £ TV 

X = 0. Hence using ( |69] ) in Lemma [3] it follows that 



CI 



tliere exists a matrix P2 > such that P2B311X = 0. 

Therefore, x'^G2 = 0, x^Csa = and B^^x = 0. Hence, 
using ( |64] | it follows that 

x'^{Gi+Gj)x^2x'^Gix 

^2X^{C2B2 + CsaBsa + G3bB3b)x = 0, 

^{Gi+Gl)x = 0, 



since, Gi + Gf > using (l72T l in Lemma |3] Therefore, 

{C2B2 + GsaBsa + Bl,Gl,)x = 0, 

MC2B2 + GsaBsa + G3aPi^Cj,)x = 0, 

in Lemma |3] Therefore, 

B2 

Bsa 



using 



[G2 Gsa] 
B2 

Bsa + P2 ^Gjf, 



^2 ^G;^ 



X = 0, 



a; = 0, 



since [G2 Gsa] is full rank using Lemma |2] Therefore, 

B2X = and (B3„ + P2^Gli,)x = 0. 

This implies that x £ M{B^a + ^2~^G^). Thus, we have 
established the second part of the lemma. Also since B2X = 



and B^bX — 0, it follows that 
if a; e TV 

Similarly, suppose that x £ TV 




0. This implies that 




, and hence 



B2 
Bsb 



X = 0. Therefore, B2X ~ 0, B^^x = and G^a; = 0. 
Hence using (l64l l in Lemma [T] it follows that 

x'^(Gi+Gf)a; = 2a;'^Gix, 

^2a:^(G2B2 + CsaB^a + C3bB3b)x = 0, 



^a;^(Gi+Gf 



0, 



since, Gi + Gf > using (iTZl i in Lemma |3] Therefore, 

x^(G2i?2 + G^bBab + B^Cl) = 0, 
^x^{G2B2 + G^bBsb + Bj,P2B3b) = 0, 
in Lemma |3] Therefore, 

X^[G2 Csb + BlP2] 

[G2 G3fa + i?3^,P2] = 0, 



using 



= 0, 



B2 
Bsb 



Hence, x £ TV 



is full rank using Lemma |2] Therefore, 

x^G2 = 0, 

^X^ [G2 G3a] = 0. 

Therefore, 




:A/' 



G2^ 
G3^. 



Hence, there exists an invertible matrix Rd such that 



'B2 




'cf 




— Rd 




Bib 





(87) 
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This completes the proof. ■ 
Lemma 6: (See |T5|) Given A e C"^" with j{A-A*)>Q 
and B e C"^" with j{B - B*) > 0, then 



det(/ - AB) ^ 0. 



(88) 



Now, we present the proof of Theorem |5] 

Proof of Theorem |5} The internal stability of the positive- 
feedback interconnection between G(s) and G{s) will be 
guaranteed by considering the closed loop system matrix 
defined in (|4]i which is given by 



A = 



A + BBC BC 
BC A 



(89) 



Here, A, B and C ai-e defined as in (I58]l-(l60ll and A, B, C, D 
are defined as in (f79]l-(l80ll. To establish internal stability, 
we show that the matrix A in (|89] | is Hurwitz; i.e., all the 
eigenvalues of A lie in the open left-half of the complex plane. 

\P - C^DC -C^C\ 



Consider T 
Lyapunov matrix, where 



P 



to be a candidate 



'Pi 








0" 


























. 








P2. 



>o, 



Pi > 0, P2 > are defined as in Lemma |3] and P > is 
defined as in Lemma |4] 

Claim 1: In the case when the matrix N in (ISTl i is a negative 
semidefinite, then T > if and only if (|83] | is satisfied. Also, in 
the case when the matrix N in dSTT l is a positive semidefinite, 
then T > if and only if and are satisfied. 

To establish this claim, we first note that since G{s) is SNI, 
it follows from Lemma |4] that P satisfies ( l86l l. This implies 
that the condition T > is equivalent to 

P - C^P»C - c^cp-^c^c > 0, 
^P - C^(P> + CP-^C'^)C > 0, 
^P - C^(P> - CA-^B)C > via (EUl in Lemmag] 

^P-C^G(0)C>0, 



'Pi 








0" 






































. 








P2. 




L'-^3f)-l 



G(0)[Ci C2 C^ia C^b] 



(91) 
(92) 



> 0. 



Furthermore, using the Schur complement of the LMI in 
it is straightforward to verify that the condition T > is 
equivalent to the conditions 



and 



-CjG(0)C2 

-gi:g(o)G2 



Pi - Gf G(0)Gi 
-G3^,G(0)Gi 



-GfG(0)G2 
-G3^,G(0)G2 

-G2^G(0)G2 
-Gf;G(0)G2 

-G2^G(0)Gi 
-Gi;G(0)Gi 



-G2^G(0)G3a' 

-G3^,G(0)G3. 



> 



(93) 



-GfG(0)G36 

P3b - G3^,G(0)G3b 
-GfG(0)G3a" 

-G3^,G(0)G3._ 

-G2^G(0)G3a' 

-G3^,G(0)G3. 

-G2^G(0)G3t" 
-G3^,G(0)G3. 



> 0. 

Furthermore 



is equivalent to 



•-^30 



where 5 is defined in 



G? 



(94) 



(95) 



(90) m 



G(0) [G2 G3a] = -S > 0, 

This is equivalent to condition 
Also, the condition (l94l i is equivalent to 



>1 




0' 

P2 




'cf 

y-ib_ 


G(0)[Gi G3b 


+ 


'cf 


G(0)[G2 G3a]S-i 


X 


''cl 


G(0) [Gi G3,.] > 0, 



>1 


0' 




"Gf 





P2 




. 3''. 



X (^G(O) - G(0) [G2 G3a] 
X [Gi G3b] > 0, 



'cl 






G(0)^ 



>1 


0' 




Gf 





P2 




. 3&. 



N [Gi Gab] > using dSB- (96) 



This condition is always satisfied in the case where N is 
negative semidefinite. Hence using ( |95] ), we can conclude that 
T > if and only if ( |83] | is satisfied in the case when N is 
negative semidefinite. 

Now in the case when N is positive semidefinite, the 
condition (|96]| can be rewritten as follows 



P/ - GjTVG/ > 0, 
Pi o' 



(97) 



where P/ = " > and G/ = [Gi Czb\. How- 

P2 

ever, using the Schur complement, WT\ is equivalent to the 
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condition 
/ 



>o, 



NiCf 
Pf 

- N^CfPJ^CjN^ > 0, 



>f 1 






. P2\ 




. 36. 



iV5 > 0, 



^ / - N^CiP^^^C^^N^ - NiC3hP2^C^bN^ > 0. (98) 
Now using ( |72] | in Lemma [3] we can define a matrix A/ as 



M + LiP^^Cf 



(99) 



so that 



M^M = C2B2 + bJ'cJ + C3B3 + BjCl 

= C2B2 + B2C2 + CsaBsa + BJ^C'^^ 

+ CsbBsb + B^,Cib- (100) 
Also using ( l68T l in Lemma [3] we can write Bi as 

= Pj^^A^Cf - L^W). (101) 
Substituting ^9} into (llOll i gives 

Bi = Pf - L^^^- - LiPf ^Cf)) 

= Pf ^Af - P^^L^M - P^^LILiP^'^CJ 
= -A^p-^Cl - P[^^L^M. (102) 

Also, from the definition of N in (ISTT l. it follows that 

[C^2 C^3a] = 0. (103) 
Therefore, Lemma |5] implies 

iv[Bj p|;]=o, 

^ iV(C2B2 + BjC2^ + CaaPsa 

+ Si;c|: + CsfaPab + i^ii^^ii)^ = 0, 
^ N{M^M)N = 0, using (fTOOl l. 

Hence, MiV = 0. Therefore, 

MN^ = 0. (104) 

Substituting this into ( |102| i implies 

PiiV^ = -AiPf ^Cf _ p-^L'^MN^, 
^BiN^ = -AiPf ^Cf7V5, 

^^TV^CiAj^^PiTV^ = -iV^CiPf ^Cfivi (105) 



Substituting (1105b into ( |98] l gives the condition 

/ + A^^CiAj-iPiA^s -iV5C36P2"VffcAf^ > 0. (106) 



This is equivalent to condition ( |84] |. Hence, in the case when 
N is positive semidefinite, it follows from (l95l l and (1106b that 
T > if and only if conditions ( |83] | and ( l84l ) are satisfied. 
This completes the proof of Claim 1. 



Now, observe that 

TA + A^T 
P ~ C^DC 
-C^C 



-c^c 
p 

T 



A + BDC BC 
BC A 



<0. 



A + BDC BC 
BC A 

P - C'^DC -C'^C 
-C^C P 

WDC + L WC 
WC L 

(107) 



Together with Claim 1, this implies that A has all its 
eigenvalues in the closed left-half of the complex plane if 
and only if conditions ( 183b and ( l84b are satisfied in the case 
when N is positive semidefinite; e.g., see Lemma 3.19 in Ii54l . 
Similarly, in the case when N is negative semidefinite A has 
all its eigenvalues in the closed left-half of the complex plane 
if and only if condition ( l83l ) is satisfied. 

In order to complete the proof of the sufficiency part of the 
theorem, we must show that if conditions ( l83b and ( |84] | are 
satisfied in the case when N is positive semidefinite, then the 
matrix A can have no eigenvalues on the joj axis. Similarly, 
we must show that if conditions ( l83b and dSSl l are satisfied in 
the case that N is negative semidefinite, then the matrix A 
can have no eigenvalues on the juj axis. 

Indeed, using Lemma|6l the fact that G(s) is NI and the fact 
that G(s) is SM, we conclude that det{I -G{ju!)G{juj)) ^ 
for all Ld > 0. This implies that A has no eigenvalues on the 
imaginary axis for lu > 0. Thus, to complete the proof, we 
will show that in the case when TV is positive semidefinite, 
conditions ( l83T l and ( l84b imply that det(^) 7^ 0. Similarly, in 
the case when N is negative semidefinite, we will show that 
conditions ( |83] | and dSSl ) imply that det(A) ^ 0. Indeed, 

det(i) = det(A) det(yl + BDC - BCA-^BC), 

= det(A)det(A + PG(0)C). (108) 

This implies that det(i) =^ if det(A + PG(0)C) ^ 0, since 
det(A) 7^ using Lemma |4] and the fact that G(s) is SNI. 
Now, define the matrix 



P2G(0)C2 P2G(0)C3a 
P3hG(0)C2 P3fcG(0)C3a 

B2 

B^h 



G(0)[C2 C3a] 



(109) 



It follows from Lemma |5] that there exists a non-singular 
matrix such that 



A-Pd 



CL 



G(0)[C2 C3a]=P<iS. (110) 



Since the matrix S is assumed to be invertible, this implies 
that the matrix A in (1 109b is invertible. 

Now, substituting (|60] | into ( 1108b . it is straightforward to 
verify that 
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det(^ + BG(O)C) 
= - det A det 

/ 





0' 








+ 







/ 


to 



X ( G(0) - G(0) [C2 C;a] A-i 
X [Ci Csb] 



\ 

= - det A det 



B2 



G(0) 



■^1 


0' 




■fir 






+ 







I 


Bat, 



G(0) - G(0) [Ga Gaa] {Rd'S.)~^Rd 



■ det A det 
' det A det 
X det + 
' det A det 



Ai 0' 

1. 
Ai 0' 
/ 

1. 
Ai 0' 
I 



X [Gi Gab 









G(0)j 



iV [Gi Gat] 



(111) 



iV [Gi Gat,] 



X det I / + I Gi Gab 

yli 0' 

' det A det 

I 



"^1 


0' 


-1 







I 




Q3a 



N 



det (/ + (GiAjT'Bi + GaftBaa) N) . 

(112) 



In the case when is positive semidefinite, ( IllOl i and ( I112l i 
imply that 



det(A + BG(0)C) 
= — det Rd det S det 

X det 
= — det Rd det S det 



^1 
/ 



I+{CiA-^B^ + C^,B^a) 
/ 



X det 



(113) 



Now using (1103b . it follows that the columns of the matrix N'^ 



are contained in the set J\f 



CI 



Hence, it follows from 



the second part of Lemma E] that {B^a + P2 ^C'sb)^^ = 0. 
This implies that 



Hence ( |113l l can be written as 

det(A + BG(O)C) 

'Ai 



— det Rd det s det 
X det 



/ 

/ + N^CiA^^BiN^ - NiC3bPi^Cjt,Ni 



(114) 



Since the matrices Rd,Ai,A are invertible and also using (|83) - 
( |84l i. dlOSI ). it follows that det (A) 7^ as required. 

In the case when N is negative semidefinite, we consider 
the matrix TV = {-N)^. Then ( 1112b impHes that 



det(^ + BG(0)C) 
= — det Rd det S det 



^1 




X det 



I - {CiA^^Bi + CsbBsa) 



det Rd det ^ det 



Ai 
/ 



X det 



I~N{CiA^^Bi+C3bB3a)N 



(115) 



Using (|103b . it follows that the columns of the matrix N are 

"cr' 



contained in the set Af 

second part of Lemma |5 
implies 



Hence, it follows from the 



that {Bsa + P2^Cj^)N = 0. This 



NCsbB^aN = ~NC3bP2'ClbN. 



(116) 



Hence, ( II 151 ) can be written as 

— det Rd det S det 



Ai 
/ 



X det 



I - NCiA-^^BiN + NC3bP2^ClbN 



(117) 



Since the matrices Rd,Ai,A are invertible and also using 
(I83ll-(l85l», ( flOSl l, it follows that dct(i) ^ as required. This 
completes the proof of the sufficiency part of the theorem. 

To complete the proof of the necessity part of the theorem, 
suppose that the positive-feedback interconnection between the 
NI transfer function matrix G{s) and the SNI transfer function 
matrix G{s) is internally stable. This implies that the matrix 
A is Hurwitz and hence has all its eigenvalues are in the open 
left-half of the complex plane. This together with Claim 1 and 
(1107b implies that conditions ( l83b and (l85b are satisfied in the 
case when N is negative semidefinite. Similarly, in the case 
when N is positive semidefinite. Claim 1 and (1107b implies 
that conditions ( [83] l and ( [84b are satisfied. This completes the 
proof of the theorem. ■ 

Now consider the following corollaries which will be used 
in proving the remaining results of the paper 

Corollary 4: Suppose that the matrix S in ( |82] | is non- 
singular and the matrix N in (ISTb satisfies N [Ci Csb] — 0. 
Also suppose that the transfer function matrix G(s), with 
the minimal state space realization (l58]l-(l59]l. is NI and the 
transfer function matrix (5(s), with the minimal state space 
realization (|79]l-(l80b, is SNI. Then the closed-loop positive- 
feedback interconnection between G{s) and G{s) is internally 
stable if and only if conditions ( l83b is satisfied. 
Proof: The proof of this corollary will proceeds in an almost 
identical fashion to the proof of Theorem |5] Indeed, we first 
state the following claim: 

Claim 2: Assume that the matrix N in (ISTb satisfies 
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N [Ci Cab] = 0, then T > if and only if ^ is satisfied. 

This claim is equivalent to Claim 1 in Theorem |5] when 
we relax the conditions on the matrix N. The proof of this 
claim is similar to the proof of Claim 1 in the proof of 
Theorem |5] since (|96] l is automatically satisfied in the case 
when N [Ci Csb] = 0. 

Also, the determinant condition in dlOSt will be automati- 
cally satisfied using the fact N [Ci Cst] = in dlllb . The 
proof of the corollary then follows as in the proof of Theorem 

m ■ 

The following corollary considers the case when n2 ~ 
and fc ^ 0; i.e., the matrix A in ( l60l l has the block diagonal 
\Ai 1 

form A = .In the case when rio = 0, the matrix 

N in (ISTl i will be given by 



(118) 



where we assume that C^G'(0)C3a is non-singular. 

We will use the following conditions in the next corollary 
which correspond to conditions (l83ll-(l85]l in Theorem |5] The 
first condition to be considered is 



(119) 



Also, for the case in which N is positive semidefinite, we will 
use the condition 



I + N^CiA^^BiN^^ - N^C3bPi^Cji,Ni > 0, (120) 

where P2 — C^a^Jbi^'ibPlib)^^ ■ Moreover, for the case in 
which is negative semidefinite, we will use the condition 



det(/ 



NCiA^^BiN + NC^bP-i^C'i^N) ^ 0, (121) 



where AT = (-A^) 2. 

Corollary 5: Suppose that the matrix C^^G{{))Cj,a is non- 
singular, fc 7^ 0, and = 0. Also, suppose that the transfer 
function matrix G(s), with the minimal state space realization 
(I58ll-(l59]l is NI and the transfer function matrix G(s), with 
the minimal state space realization in (l79]l-(l80ll, is SNI. If 
A^ in (II I8I 1 is positive semidefinite, then the closed-loop 
positive-feedback interconnection between G{s) and G{s) is 
internally stable if and only if conditions (11 19t and ( I120l l are 
satisfied. Also, if A^ in (1118b is negative semidefinite, then the 
closed-loop positive-feedback interconnection between G{s) 
and G{s) is internally stable if and only if conditions ( 11 \9\ 
and ( 1121b are satisfied. 

Proof: Note that for the system (l58Tl-(l60]l corresponding to 
the case of this corollary, the conditions (|83Tl-(l85Tl in Theorem 
|5] reduce to conditions ( 11 19b - (ll21b . Then the proof of the 
corollary proceeds in an identical fashion to the proof of 
Theorem |5] for the special case being considered, where the 
matrix P defined in ( |90l l becomes a matrix of the form 
"Fi ^ 



P 







^0 P2J 

The next corollary considers the case when n 7^ and k — 
0; i.e., the A matrix in the minimal state realizatioriof G{s) 
, _ Ai 

b-(l60b has the block diagonal form A— a 

A'} 



case, when n ^ {) and fc = 0, the matrix A^ in (ISTb will be 
given by 



N^V{G{0),C2) 



(122) 



where the matrix Cj(5(0)C2 is assumed to be non-singular. 

We will use the following conditions in the next corollary 
which corresponds to conditions (ISJb-dSSTl in Theorem |5] The 
first condition to be considered is 



C2^G(0)C2 < 0. 



(123) 



Also, for the case in which A^ in (1122b is positive semidefinite, 
we will use the condition 



I + N^^GiA'^^Bim > 0. 



(124) 



Moreover, for the case in which A^ in (1122b is negative 
semidefinite, we will use the condition 



det(/ - NCiA^^BiN) ^ 0, 



(125) 



where A^ = i-N)^. 

Corollary 6: Suppose that CjG(0)C2 is non-singular, 
n2 ^ 0, and fc = 0. Also, suppose that the transfer function 
matrix G(s), with the minimal state space realization (|58]|-(|59]| 
is NI and the transfer function matrix G(s), with the minimal 
state space realization in (|79]l-(l80ll, is SNI. If A^ in ( 1122b is 
positive semidefinite, then the closed-loop positive-feedback 
interconnection between G(s) and G(s) is internally stable 
if and only if conditions (1123b and (1124b are satisfied. Also, 
if A^ in (1122b is negative semidefinite, then the closed-loop 
positive-feedback interconnection between G(s) and G(s) is 
internally stable if and only if conditions (1123b and (1125b are 
satisfied. 

Proof: First note that for the system (I58]l-(l60ll corresponding 
to the case of this corollary, the conditions (l83]l-(l85ll in 
Theorem |5] reduce to conditions (I123t - (ll25b . Then the proof 
of the corollary proceeds in an identical fashion to the proof 
of Theorem |5] for the special case being considered, where 
the matrix P defined in ( |90b becomes a matrix of the form 
>i 0" 



P 







Now to present the proof of the main results of the paper, we 
will use Theorem|5] Corollaries HMJand the following lemmas. 
We first show that any NI system can be transformed to the 
block diagonal form given in (l58]l-(l60ll. 

Lemma 7: Any NI system with transfer function matrix 
G(s) and minimal state space realization ( fSSl l. ( |59] l, can be 
transformed to the block diagonal form given in (l60t . 

Proof: Suppose that the transfer function matrix G(s), 



with a minimal state space realization 



A 


B 


G 






NI. 



It follows from Theorem 2.1.1 in |61| that we can find a 
non-singular state space transformation matrix T such that 
the matrix T~-^AT is in Jordon block diagonal form and 

■ T-^B ^ 



the realization 



T-^AT 



CT 







is minimal. Also, we 



can choose this transformation so that the Jordon blocks 
of T^^AT are ordered according to the magnitudes of the 
corresponding eigenvalues of the matrix A, such that the 
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last blocks correspond to the zero eigenvalues of A if they 
exist. Furthermore, this transformation can be chosen so that 
the Jordan blocks corresponding to the zero eigenvalues are 
ordered according to increasing order of the Jordan blocks. 
Also, a further transformation can be applied so that the matrix 
A3 corresponding to the Jordan blocks of order two is of the 
form ( leTT i. 

Now, we claim that if G(s) is NI, then there are no Jordan 
blocks corresponding to zero eigenvalues of order greater than 
or equal to three. To prove this claim, suppose that there is a 
Jordon block of A corresponding to a zero eigenvalue of order 
greater than or equal to three. This together with the minimal- 
ity of the realization implies that G3 = lim s^G{s) ^ 

s >0 

which contradicts the NI definition. Thus the zero eigenvalues 
of A will only have Jordon blocks of order one or two. From 
this, it now follows that the matrix T^^AT will be of the form 
( |60] |. This completes the proof of the lemma. ■ 

The next lemma is a technical lemma, which will be used 
in order to prove our results. 

Lemma 8: For any full rank matrices A, B, C and D which 
satisfy AB = CD where A e R"><^S e W''\C £ 
]^nxr^ 13 e M''^" and n > r, there exists an invertible matrix 
R such that A = CR and B = R-^D. 

Proof: Since B is of full rank, and n > r, AB — CD 
implies 

ABB^ = CDB^, 
^A = CDB^{BB^)-\ 
^A^CR 

where R ^ DB^{BB^)-^. 

To show that R is nonsingular, suppose that R is singular. 
Then there exists a non-zero nxl vector x such that Rx ~ 0. 
This implies that Ax = which contradicts the fact that A 
is a full rank. Hence, that there exists a nonsingular matrix R 
such that A = CR. Also, since C is of full rank and n > r, it 
follows that C has a left inverse, which implies that RB = D; 
i.e., B — R^^D. This complete the proof. ■ 

We are now ready to give the proof of Theorem [T] 

Proof of Theorem [7} Lemma [T] shows that any NI system 
can be represented in the block diagonal form (l5Ft -(l6Qt. This 
implies that we only need to show the equivalence of the 
assumptions and the conditions (l83Tl-(l85Tl in Theorem |5] and 
the assumptions and the conditions (fT6ll-(fT8l) in this theorem. 

First, it is straightforward to verify that the condition k ^ {] 
is equivalent to the condition G2 ^ 0. Also, it follows from 
(l63T l and (|9]l that there exists a full rank matrix J such that 



CsaB^b — JJ 



(126) 



Also, it follows from Lemma |2] and Lemma |8] that there 
exists an invertible matrix X such that C^a — JX and 
B^b = X^^JT. We let P2 = C^^B^^iB^bBj^)-^ and note 
that i?3f,i3j{, is invertible since B^b is of full rank. Then 
Lemma [3] implies that P2 is symmetric and also we obtain 



(127) 



In the case when N is positive semidefinite, the definition of N 
implies that N [C2 Csa] = 0, and hence [C2 Csa] = 



+ 



0. Using (|64ll, it follows that Gi = [G2 Csa] 
CsbB^b, which implies 

N^Gi^N^C2 Gsa] 

^ iV^Gi = N^CsbBsb since [C2 Csa] = 0, 

N^-GiJ = N^^CzbX-^J^J, 
^ N^GiJiJ^J)-^ = N^C3bX-\ (128) 



N^CsbBsb, 



Substituting (|65j, ( |127t and ( |128t into condition in 
Theorem |5] it follows that this condition can be rewritten as 

/- (^mCoN^ + mCiJiJ^ J)-^{J^ J)-'^ J'^GJ Ni^ 

> 0. (129) 

Similarly, in the case when N is negative semidefinite, it 
follows that N [C2 Csa] = 0. This implies that 



NGi^N[C2 Csa] 



+ NCsbBsb, 
Gaa] = 



B2 

Bsa 

NGi = NC^bBzb since N [C2 Gga] = 0, 
^ NGiJ = NC3bX~^j'^J, 

=> NGiJ{J^jy^ = NC3bX~\ (130) 

Substituting (|65]l, ( |127t and dlBOt into condition (|85]l in 
Theorem |5] it follows that this condition can be rewritten as 

det (/ + (^NGaN + NGiJ{J^ J)-\j'^ j'^Cj N^^ 

^0. (131) 

Now, using Lemma [T] and substituting for Gi and G2 from 
dill, (It) in the Hankel matrix defined in ( fTOl l. it follows that 



where 



r = 

~ U2 
A = 

I ^ 

G= [G2 



Gi 


G2' 




" C ' 


G2 







CA 



[B AB] (132) 





A3 
G3I 





"0 


' 




'B2 






A3_ 


; B^ 











B3 



Using the SVD in ([IB and (fT32l ). it follows that 

C 
CA 



[13 AI3] 



Ui 
U2 



(133) 



Using Lemma |2] and Lemma |8] it follows that there exists a 
nonsingular matrix R such that 



(134) 









' c ' 




" G " 


u = 


U2 




CA 


R = 


CA 



where 

C^CR, 
A = R-^AR. 

This implies that A^ = R-^A^R = since A^ = 0. 
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It follows that UA 
implies 



Using ([T2]i and (|135t . it follows that 

AfiA) = span{V2}- 
"0 

A-, 

Also, since A 



CA 








U2 






n 








"c/2' 




= [/ 







which 



(135) 





A3 



(136) 
, it follows that 



0/ 


7 0' 

Af{A) = span{ 1}. (137) 


Now observe that we can write the matrix [C2 Csa] as 

7 0' 

[C2 Caa] - C / . (138) 


Also, observe that Ax = if and only if 

R-^ARx = 
^ARx = 
^Rx eJV{A). 



Hence, J\f{A) ^ RM{A). Therefore it follows from (IT36^ and 
(fT37] i that 



i? span{V2} = spari{_RV2} = span{ 



I 
/ 




This implies that there exists a nonsingular matrix R such that 

7 0' 

/ ^RV2R. (139) 




Substituting (|139l l into (|138t and using (|134t impUes 

[C2 Csa] = Ci?T/2i? = CV2R - {/iT^2/? = /^i?, (140) 

where F = U1V2 as in (fTST l. Substituting (|140| i into the matrix 
(ISTT i and (|82]i implies that 

N = (g(0) - G{0)FR{R^F'^G{0)FR)-^It^F'^G{0)^ 

= G(0) - G{0)F {F'^G{0)Fy'^ F'^G{0) 

= Nf, (141) 



where Nj is defined as into (O. Substituting ( 1141b in (|129l l 
and ( II3II 1 implies that conditions ( |84] | and dSSl l in Theorem 
|5] are equivalent to conditions ([TtT i and ( fTSl l in the theorem 
respectively. 



Also, ( |140| i implies that 

G(0) [G2 G3a] = lfF^G(Q)FR. 



•-^30 



It follows that condition ( [83T l in Theorem |5] is equivalent 
to condition ( fT6] l in the theorem since i? is invertible. This 
completes the proof of the theorem. ■ 
Proof of Corollary [7} In order to prove this corollary, we 
show that the stability conditions and the assumptions in 
Corollary |5] are equivalent to the stability conditions and the 
assumptions in this corollary. First, it is straightforward to 
verify that the conditions fc and ri2 ~ are equivalent 
to the conditions G2 ^ and Gi = 0. Also, using ( |63] | and 
the decomposition in (|9]l, it follows that GsaB^b = JJ^, and 
hence Lemma [8] implies that there exist an invertible matrix X 
such that Gsa = JX. This implies that the matrix A'^2 in (O 
is equal to the matrix N in (11 18b . Also, since C^a = JX and 
X is invertible, it follows that condition ( |20] | in Corollary [T] 
is equivalent to condition (11 191 l in Corollary |5] Since Gi = 0, 
it follows that 



NGi - NiCsaB^a + CsbBsb) = 0, 
^ NC3bB3b = 0, since NGsa = 0, 
NCsb — 0, since B^b is of full rank. 



(142) 



This implies that N^C^b = in the case when N is positive 
semidefinite. Using the fact that Gq = —CiA^^Bi from 
Lemma [T] it follows that condition ( |2TI ) in Corollary [T] is 
equivalent to condition ( 11201 ) in Corollary |5] Also, in the case 
when N is negative semidefinite ( 1142b implies that NC^b = 0. 
Using the fact that Gq — —CiA^Bi from Lemma [T] it 
follows that condition (l22l i in Corollary [T] is equivalent to 
condition (1121b in Corollary |5] This completes the proof of 
the corollary. ■ 
Proof of Theorem ^ In order to prove this theorem, we 
first show that N{G2) ^ N{G'^) implies the condition 
N [Gi Gab] = in Corollary g] in the case when Gi — 0. 
Indeed, suppose that M{G2) Q N{G'^). This implies that 



7^(G2) 3 7^(Go) 



(143) 



where 7^(-) denotes the range space of a matrix. 

Since G2 = JJ^ and J is of full rank, it follows that 



7^(JJ^) D 7^(Go) 



(144) 



which implies that there exist a matrix Q such that Go ~ JQ- 
Then, we consider the matrix N defined as 

N = G(0) - G(0)G3a(G3^aG(0)G3a)-V3'lG(0), (145) 

which is the form of the matrix N in Corollary |4] in the case 
in which Gi = 0. This implies that 

iVGo = NJQ = 0, since NJ = 0, 

and hence from Lemma [T] it follows that 

NCiA-^Bi ^ 0. (146) 

Using a similar calculation as in equation ( 1102b in the proof 
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of Theorem |5] (|146t implies that 

NCiA]:\AiP^^C^ - P{^l1m)N ^ 0, 
^ A^CiPf ^Cf = using (fT04l l. 
^ iVCi = 0. (147) 

Also, since Gi = CsiJs = CsaB^a + CzhB-^b = 0, it follows 
that 

NC^aB^ia + ^^36^36 = 0, 

=^ NCsbBsb = 0, since ArC3a = 0, 

iVCsb = since B36 is of full rank. (148) 

Using (fl47] i and (fT48] i. it follows that N [Ci C^b] = 0. This 
implies the assumptions in Corollary |4] are satisfied in the 
case when Gi = 0. Also, as in the proof of Theorem [T] the 
condition ( l20l i reduces to condition ( |83] | in Corollary |4] ■ 

Proof of Corollary |2} Since G2 is assumed to be positive 
definite in this corollary, it follows that the matrix J in 
(|9]l is invertible. Then, the condition dZSl l reduces to the 
condition (5(0) < and the corollary follows immediately 
from Corollary |2] ■ 

Proof of Theorem \3}i In order to prove this theorem, we 
show that the stability conditions and the assumptions in this 
theorem are equivalent to the stability conditions and the 
assumptions in Corollary |6] First, it is straightforward to verify 
that the conditions ri2 7^ and fc = are equivalent to the 
conditions Gi ^ and G2 — 0. Using Lemma |8] and the fact 
that Gi = G2i?2 from Lemma [T] it follows that there exists 
an invertible matrix R such that G2 = FiR, where the matrix 
Fi is given in ( |23] |. This implies that the matrix A^i in (l24l i 
is equal to the matrix N in (I122l l. Also, since G2 = FiR and 
R is invertible, it follows that condition dZST l in this theorem 
is equivalent to condition ( 1123b in Corollary |6] Finally, using 
the fact that Go — —CiA^^Bi from Lemma [T] it follows 
that conditions ( |26] | and dZTl i in this theorem are equivalent to 
conditions (1124b and ( 1125b in Corollary |6] respectively. This 
completes the proof of the theorem. ■ 



Proof of Theorem |5} In order to prove this theorem, first 
we show that 7V(Gf) C 7V(Gj5^) implies the condition 
N [Gi Gab] = in Corollary H in the case when G2 = 0. 
Indeed, suppose that 7V(Gf ) C M{Gq). This implies that 

7^(Gl) D 7^(Go). (149) 

Since Gi ~ C2B2 from Lemma [T] and B2 is of full rank using 
Lemma [21 it follows that TZ{C2) = Tl{Gi). Using ( fT49] l. it 
follows that 

niC2) 2 nGo), (150) 

which implies that there exist a matrix Q such that Go = G2Q. 
Then, we consider the matrix N defined as 

N = G(0) - (5(0)G2(G2^G(0)G2)~^G2^G(0), (151) 

which is the form of the matrix N in Corollary |4] for the case 
in which G2 = 0. This implies that 

NGo = NC2Q = 0, since iVG2 = 

and hence from Lemma [T] it follows that 

NCiA^^Bi ^ 0. (152) 

Using a similar calculation as in equation ( 1102b in the proof 
of Theorem |5] (1152b implies 

NCiA^\AiP^^Gl - P^^lIM)N ^ 0, 
^ NCiP^^ClN = 0, using 
^ NCi = 0. 

This implies the assumptions in Corollary |4] are satisfied in 
the case when G2 = 0. Also, as in the proof of Theorem [T] 
condition (|25] | reduces to condition ( |83] | in Corollary |4] ■ 
Proof of Corollary\3^ Since Gi is assumed to be invertible in 
this corollary, it follows that the matrix Fi in ( |23] | is invertible. 
Then, the condition (l25b reduces to the condition (5(0) < 
and the corollary follows immediately from Corollary 5] B 



